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ABSTRACT 

The  Army  is  currently  developing  acoustic  sensor  systems  that  will  provide  extended  range  surveillance, 
detection,  and  identification  for  force  protection  and  tactical  security.  A  network  of  such  sensors  remotely 
deployed  in  conjunction  with  a  central  processing  node  (or  gateway)  will  provide  early  warning  and 
assessment  of  enemy  threats,  near  real-time  situational  awareness  to  commanders,  and  may  reduce 
potential  hazards  to  the  soldier.  In  contrast,  the  current  detection  of  chemical/biological  (CB)  agents 
expelled  into  a  battlefield  environment  is  limited  to  the  response  of  chemical  sensors  that  must  be  located 
within  close  proximity  to  the  CB  agent.  Since  chemical  sensors  detect  hazardous  agents  through  contact, 
the  sensor  range  to  an  airburst  is  the  key-limiting  factor  in  identifying  a  potential  CB  weapon  attack.  The 
associated  sensor  reporting  latencies  must  be  minimized  to  give  sufficient  preparation  time  to  field 
commanders,  who  must  assess  if  an  attack  is  about  to  occur,  has  occurred,  or  if  occurred,  the  type  of 
agent  that  soldiers  might  be  exposed  to.  The  long-range  propagation  of  acoustic  blast  waves  from  heavy 
artillery  blasts,  which  are  typical  in  a  battlefield  environment,  introduces  a  feature  for  using  acoustics  and 
other  sensor  suite  technologies  for  the  early  detection  and  identification  of  CB  threats.  Employing 
disparate  sensor  technologies  implies  that  warning  of  a  potential  CB  attack  can  be  provided  to  the  solider 
more  rapidly  and  from  a  safer  distance  when  compared  to  current  conventional  methods.  Distinct 
characteristics  arise  within  the  different  airburst  signatures  because  High  Explosive  (HE)  warheads 
emphasize  concussive  and  shrapnel  effects,  while  chemical/biological  warheads  are  designed  to  disperse 
their  contents  over  immense  areas,  therefore  utilizing  a  slower  burning,  less  intensive  explosion  to  mix 
and  distribute  their  contents.  Highly  reliable  discrimination  (100%)  has  been  demonstrated  at  the 
Portable  Area  Warning  Sun’eillance  System  (PAWSS)  Limited  Objective  Experiment  (LOE)  conducted  by 
Joint  Project  Manager  for  Nuclear  Biological  Contamination  Avoidance  (JPMNBC  CA)  and  a  matrixed 
team  from  Edgewood  Chemical  and  Biological  Center  (ECBC)  at  ranges  exceeding  3km.  The  details  of 
the  field-test  experiment  and  real-time  implementation/integration  of  the  standalone  acoustic  sensor 
system  are  discussed  herein. 


1.  INTRODUCTION 

The  Army  is  currently  developing  acoustic  sensor  systems  that  will  provide  extended  range  surveillance, 
detection,  and  identification  for  force  protection  and  tactical  security.  A  network  of  such  sensors  remotely 
deployed  in  conjunction  with  a  central  processing  node  (or  gateway)  will  provide  early  warning  and 
assessment  of  enemy  threats,  near  real-time  situational  awareness  to  commanders,  and  may  reduce 
potential  hazards  to  the  soldier.  In  contrast,  the  current  detection  of  CB  agents  expelled  into  a  battlefield 
environment  is  limited  to  the  response  of  chemical  sensors  that  must  be  located  within  close  proximity  to 
the  CB  agent.  Since  chemical  sensors  detect  hazardous  agents  through  contact,  the  sensor  range  to  an 
airburst  is  the  key-limiting  factor  in  identifying  a  potential  CB  weapon  attack.  The  associated  sensor 
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reporting  latencies  must  be  minimized  to  give  sufficient  preparation  time  to  field  commanders,  who  must 
assess  several  situations  of  the  enemies’  strategy,  as  well  as  factoring  the  exposure  time  and  types  of 
agents  the  soldiers  could  encounter  on  the  battlefield.  The  long-range  propagation  of  acoustic  blast  waves 
from  heavy  artillery  blasts,  which  are  typical  in  a  battlefield  environment,  introduces  a  feature  for  using 
acoustics  and  other  disparate  sensor  technologies  for  the  early  detection  and  identification  of  CB  threats. 
Employing  disparate  sensor  technologies  implies  that  warning  of  a  potential  CB  attack  can  be  provided  to 
the  solider  more  rapidly  and  from  a  safer  distance  in  comparison  to  conventional  methods  lessening 
exposure  time  for  the  soldiers.  This  capability  facilitates  the  need  for  classifying  the  types  of  rounds  that 
have  burst  in  a  specified  region  in  order  to  give  both  warning  and  provide  identification  of  CB  agents 
found  in  the  area.  The  additional  information  incorporated  into  the  common  operating  picture  will  help 
field  commanders  (more)  accurately  plan  and  execute  measures  that  must  be  taken  in  order  to  ensure  the 
proper  precautions  to  a  mission  [7].  Such  a  warning  system  would  also  help  reduce  the  time-consuming, 
manpower  intensive  and  dangerous  tasks  associated  with  identifying  the  airburst,  allowing  a  field 
commander  to  make  rapid  and  accurate  judgments  that  insure  greater  safety  and  lessen  exposure  for  the 
soldiers. 

This  paper  supports  acoustic  sensor  technologies  to  differentiate  among  conventional,  i.e.  HE,  and  CB 
artillery  blasts  using  robust  characteristics  that  remain  unchanged  with  long-range  wave  propagation  and 
degradation,  firing  time  and  detonation  (either  air  or  ground).  The  outcome  illustrated  herein  is  based  on 
an  analysis  conducted  using  acoustic  signature  data  collected  during  two  independent  data  collection 
exercises  and  a  third  exercise  that  was  conducted  to  demonstrate  the  algorithms  in  a  real-time 
environment.  The  principal  set  of  artillery  blast  data  was  collected  at  YUMA  Proving  Ground  on  11-12 
August  1995  in  a  field  test  experiment  conducted  by  the  National  Center  of  Physical  Acoustics  (NCPA)  at 
the  University  of  Mississippi  in  collaboration  with  Acoustic  Center  of  Excellence  (ACOE)  at  Picatinny 
Arsenal.  The  reason  for  the  experiment  was  to  study  if  CB  artillery  blasts  could  be  distinguished  from  HE 
blasts  via  acoustic  and  seismic  sensor  technologies.  The  field  test  consisted  of  point  detonation  and/or 
airburst  of  three  categories  of  explosive  rounds;  two  categories  of  artillery  rounds  considered  analogous  to 
chemical  and  biological  warheads,  and  generic  HE  rounds  most  commonly  used  in  the  battlefield.  The 
data  set  comprised  of  8  ground  impacts  from  HE  rounds,  23  ground  impacts  from  a  Type  A  CB  round  and 
8  air-bursts  from  a  Type  B  CB  round.  The  data  was  collected  at  seven  sensors  sites  as  shown  in  Figure  1. 
The  rounds  were  fired  from  a  self-propelled  howitzer  situated  several  kilometers  to  the  west  of  a  linear 
array  of  5  sensor  locations  designated  along  an  axis  offset  from  the  firing  radial  by  700m  (cf.  Figure  1). 
Sensor  site  S3  coincides  with  the  ordinate  axis  of  the  estimated  projectile  impact  point  at  5km  from  the 
firing  position.  Microphones  were  emplaced  at  sensor  sites  S2,  S4,  SI  and  S5  at  a  distance  of  500m, 
600m,  1000m  and  1200m  from  S3,  respectively.  Geophones  accompanied  the  microphones  at  positions  S3 
and  S5,  and  a  flash  detector  was  emplaced  at  position  S6,  2000  meters  from  S3.  An  acoustic  transducer 
was  positioned  at  site  S7,  400m  from  S3  and  in  line  with  the  estimated  projectile  impact  point. 


Figure  1  -  Sensor  Deployment  at  Yuma  Proving  Grounds 


A  preliminary  investigation  conducted  by  researchers  at  the  National  Center  for  Physical  Acoustics  has 
shown  distinctive  differences  between  airburst  signatures  from  the  three  types  of  rounds,  and  a  subsequent 
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study  identifies  potential  features  within  the  blast  waves  that  could  be  used  to  discriminate  between  HE 
and  CB  blasts  acoustically  [1],  In  this  study,  the  rapid  ignition  of  the  HE  rounds  was  found  to  generate 
larger  amplitude,  shorter  duration  signatures  in  contrast  to  the  longer  duration  signals  representative  of  the 
slower  burning  phosphorus,  ignited  in  the  Type  A  shell  upon  ground  impact.  It  was  determined  that 
accurate  range  information  from  the  multiple  sensor  positions  to  the  blast  event  location  would  permit  the 
use  of  amplitude  alone  to  identify  a  CB  round  from  its  conventional  HE  counterpart.  Discriminating 
attributes  such  as  high  frequency  precursor  signals  produced  by  the  passing  of  supersonic  shrapnel 
elements  from  the  exterior  of  HE  rounds  were  identified  prior  to  the  main  blast  wave  and  noticeable  at 
ranges  up  to  1  km,  absent  from  CB  rounds. 

A  second  set  of  experimental  data  was  collected  at  Dugway  Proving  Ground  (DPG)  to  examine  if  acoustic, 
seismic,  and  imagery  sensors  could  be  exploited  in  order  to  distinguish  conventional  HE  projectiles  from 
those  filled  with  either  chemical  and/or  biological  agents  [2].  During  this  test  exercise,  two  hundred  and 
sixty  rounds  were  fired,  approximately  half  of  which  were  fused  to  detonate  on  impact,  while  the  other 
half  were  fused  to  detonate  in  the  air.  The  HE  rounds  contained  a  larger  amount  of  explosives  when 
compared  to  the  simulated  chemical  rounds  specifically  developed  to  detonate  and  release  a  fill  consisting 
of  polyethylene  glycol  and  water.  The  data  examined  in  this  paper  was  collected  from  three  sensor  sites 
arranged  in  a  triangle  and  numbered  according  to  range  from  the  detonation  zone.  Site  1  and  site  3  were 
situated  on  the  same  radial  line  from  the  expected  detonation  area  at  ranges  of  540m  and  1665m 
respectively.  Site  2  was  positioned  1026m  radially  outward  and  5 14m  to  the  left  of  said  location. 

Numerous  studies  have  been  conducted  that  have  created  a  wealth  of  research  concerning  the 
discrimination  of  HE  and  CB  rounds  using  disparate  sensor  technologies  [2].  Both  time  and  frequency 
domain  approaches  were  isolated  to  obtain  features  useful  in  the  discrimination  between  HE  and  CB 
rounds.  Time  domain  discriminating  features  were  identified  within  the  blast  envelope,  including  the 
shrapnel  noise  reported  in  [1],  and  frequency  domain  features  such  as  the  instantaneous  and  maximum 
bandwidth. 

A  third  field  exercise/data  collection  was  conducted  to  test  and  benchmark  the  real-time  performance  of 
the  CBRN  discrimination  algorithm  developed  by  ARDEC.  The  JPM  NBC  CA  office  provided  ARDEC 
the  opportunity  to  use  the  PAWSS  system  as  a  test  platform  for  the  algorithm  described  herein  developed 
by  the  ACOE  [6].  The  PAWSS  platform  contained  MET  sensors,  Mid  Range  Radar,  Stand-OFF  CB 
sensor,  Long  Range  Radar,  Point  CB  sensors,  Battlefield  Warning  software  and  Surveillance  Systems 
providing  a  one  stop  CBRN  detection  solution.  The  artillery  airburst  portion  of  the  experiment  was  held 
between  20-22  June  2005  as  part  of  the  PAWSS  LOE  conducted  by  JPM  NBC  CA  and  their  team 
matrixed  from  ECBC. 

The  most  distinctive  characteristics  of  each  blast  type  are  easily  identified  in  Figure  2;  the  high  frequency 
precursors  from  the  HE  blast  prior  to  the  main  blasts  as  shown  in  Figure  2a,  the  comparatively  slow  bum 
and  weak  under  pressure  associated  with  the  Type  A  CB  blast  in  Figure  2b,  and  a  short  duration  pulse  that 
precedes  the  slow  bum  of  agents  within  the  Type  B  Base  ejection  rounds  shown  in  Figure  2c.  We  show 
that  range  invariant  features  found  within  the  multiresolution  decomposition  of  the  corresponding 
signatures  in  an  attempt  to  further  promote  disparate  sensor  technology  for  CB/HE  discrimination.  The 
subsequent  growth  of  a  feature  space  representation  for  classification  is  based  upon  data  acquired  during 
separate  data  collection  exercises,  that  when  combined,  introduce  variations  within  the  intraclass  and 
interclass  attributes  for  the  HE  and  Type  A  CB  artillery  blasts.  This  is  principally  due  to  the  different 
types  of  payloads  used  to  simulate  CB  rounds,  differences  in  the  local  environment  and  data  acquisition 
hardware  used,  and  variations  in  range  and  aspect  angle  of  the  blasts  with  respect  to  sensor  emplacement 
positions.  We  show  that  information  extracted  via  different  levels  of  decomposition  produce  repeatable 
properties  that  can  be  combined  to  ensure  reliable  discrimination  between  HE  and  CB  blasts. 
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Tima  loocond] 


Figure  2  -  Acoustic  Signatures  from  (a)  HE  blast  recorded  at  DPG,  (b)  Type  A  Simulated  CB 
Round  burst  recorded  at  DPG  and  (c)  Air  burst  for  the  Type  B  White 
Phosphorus  Base  Injection  round  recorded  atYPG 


2.  WAVELET  FEATURE  EXTRACTION  FOR  CLASSIFICATION 

Wavelet  analysis  is  used  to  identify  distinct  feature  sets  that  remain  consistent  for  a  given  class  of  artillery 
round,  and  do  not  degrade  significantly  with  long-range  propagation.  Figure  2  shows  the  complex,  non¬ 
stationary  response  of  the  signatures  to  be  evaluated,  which  are  categorically  poor  candidates  for  feature 
extraction  and  segmentation  via  Fourier  analysis  or  the  short-time  Fourier  Transform.  On  the  other  hand, 
the  non-stationary,  transient  and  often  oscillatory  nature  of  these  signals  is  readily  characterized  with 
wavelet  bases  that  effectively  capture  the  time-frequency  distribution  of  such  signal  components. 
Wavelets  are  better  suited  for  analyzing  transient  signals,  as  they  are  well  localized  in  time  and  in 
frequency,  and  are  able  to  take  into  account  the  scale  of  multiple  signal  components  [3].  To  this  end,  the 
wavelet  transform  can  provide  a  scalable  time-frequency  representation  of  artillery  blast  signatures  and 
provides  details  that  are  not  readily  found  using  conventional  signal  processing  techniques.  Noise 
components  that  are  introduced  into  the  blasts  waves  prior  to  being  received  at  the  sensor  are  also  a 
serious  impediment  to  isolating  features  that  can  be  used  to  efficiently  characterize  the  signal. 
Multiresolution  analysis  will  also  isolate  this  noise  and  remove  it  from  the  baseband  signal  of  interest;  a 
process  referred  to  as  wavelet  denoising.  Wavelet  denoising  is  optimal  in  that  noise  components  are 
removed  from  signal  components  regardless  of  the  frequency  content  of  the  signal,  which  is  by  far  more 
effective  than  conventional  filtering  methods. 

2.1  Discrete  Wavelet  Transform  and  Multiresolution  Analysis 

The  discrete  wavelet  transform  (DWT)  is  derived  from  subband  filters  and  based  on  a  multiresolution 
decomposition  of  a  signal  to  give  a  coarser  and  coarser  approximation  to  the  original  signal  by  removing 
high  frequency  detail  at  each  level  of  decomposition  [4].  In  other  words,  the  wavelet  transform  is  a 
multiresolution  transform  that  maps  low  frequency  information  of  signals  into  a  coarsely  sampled 
subspace  and  maps  high  frequency  information  into  more  finely  sampled  subspaces.  The  DWT  is  defined 
by  the  scaling  function 
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and  a  wavelet  function: 


,L~  1 


4*) = 2/2  X  _  *) 


k=o 


yL-l 

w{x)  =  2/2  ^  gk+A2x  ~  k) 


k=0 


(1) 


(2) 


where  hk  and  gk  are  analysis  filters.  It  is  important  to  choose  the  appropriate  filter,  in  order  to  retain  the 
characteristics  of  the  transient  signals.  .  The  choice  of  the  wavelet  scaling  function  should  have  properties 
similar  to  the  original  signal  because  the  quality  of  wavelet  decomposition  depends  significantly  on  the 
ability  to  approximate  the  signal  with  wavelets  [3].  The  output  of  the  wavelet  transform  shows  the 
correlation  between  the  signal  and  wavelet  as  a  function  of  time.  The  easiest  method  for  selecting  a 
wavelet  is  to  copy  a  signal’s  time-frequency  behavior,  however,  most  transient  signals  cannot  be  used  as 
wavelet  basis  functions  since  they  resemble  exponentially  damped  sinusoids  and  do  not  possess  a  zero 
mean.  The  wavelet  basis  used  in  this  paper  is  the  db5  wavelet  defined  by  Daubechies  [3]  that  has  the 
scaling  function  and  translation  function  shown  in  Figure  3(a),  and  3(b)  respectively. 


Scaling  Function  phi  Translation  Function  psi 

1.2 - , - * - , - , - , - , - , - , -  1.5  - - , - t - , - , - , - , - 


Figure  3  -  Daubechies  wavelet  n  =5  (a)  Scaling  Function  (b)  Translation  Function 

Note  that  scaling  function  strongly  resembles  the  blast  signature  for  the  HE  and  "Type  A  CB  signals" 
shown  in  Figure  3a  and  3b,  as  well  as  the  short  pulse  preceding  the  main  blast  in  Figure  2c.  A 
multiresolution  filter  bank  using  low  pass  and  high  pass  wavelet  filters  as  shown  in  figure  4  is  used  to 
implement  the  DWT  and  to  decompose  an  input  signal  into  different  frequency  bands  is  used  to  implement 
the  DWT.  At  each  level  of  decomposition,  the  high  pass  filter  defined  in  (2)  produces  the  detail 
information,/),, ,  while  the  low  pass  filter  associated  with  the  scaling  function  in  (1)  produces  the  coarse 
approximations  An  .  Often  referred  to  as  the  Mallat  algorithm  a  process  of  successive  lowpass  and  high 
pass  filtering  of  the  input  signal  is  used  to  implement  the  DWT  [5].  The  resulting  hanks  of  dyadic 
multirate  filters  split  the  input  signal  into  subbands  at  each  level  of  decomposition  produces  a  signals  with 
a  subset  of  frequencies  spanning  half  of  the  original  frequency  band.  If  the  original  signal  is  sampled  at  a 
frequency  of  fs  Hz,  then  the  output  of  the  first  high  pass  filter  which  is  the  first  detail  coefficient  Dx , 
captures  the  band  of  frequencies  between  fs/2  and  fs/4.  In  the  same  manner,  the  highpass  filter  in  the 
second  stage  captures  signal  components  with  a  bandwidth  between  fs/4  and  fs/8,  and  so  on.  In  this  way, 
arbitrary  resolution  in  frequency  can  be  obtained.  Since  the  input  signal  at  each  stage  of  decomposition 
contains  the  highest  frequency  equivalent  to  twice  that  of  the  output  stages,  it  can  be  sampled  at  half  the 
original  sampling  frequency,  thus  discarding  half  the  samples  with  no  loss  of  information.  This 
decimation  by  a  factor  of  2,  reduces  the  time  resolution  of  the  entire  signal  as  the  input  signal  is 
represented  by  half  of  the  total  number  of  samples  and  effectively  doubles  the  scale. 


RTO-MP-SET-107 


8-5 


UNCLASSIFIED/UNLIMITED 


UNCLASSIFIED/UNLIMITED 

Reliable  Classification  of  High  Explosive  and 
Chemical/Biological  Artillery  Using  Acoustic  Sensors 


ORGANIZATION 


Figure  4  -  A  Multirate  Filter  Bank  used  as  a  Five  Level  Wavelet  Decomposition  Tree 


2.2  Wavelet  Features  for  Classification 

Figures  5  (a-c)  show  blast  waves  from  the  three  unique  artillery  rounds  of  interest,  followed  by  four 
signals  derived  from  the  multiresolution  decomposition  with  the  coefficients  at  level  5  and  detail  signals  at 
level  5,  4,  and  3  respectively.  The  proposed  feature  space  is  comprised  of  primitives  derived  from  the 
normalized  energy  distributions  within  the  details  at  level  5,  4,  and  3,  centered  about  the  maximum  value 
of  the  blast  wave.  A  fourth  feature  is  obtained  from  the  coefficients  at  level  5.  Features  such  as  the  rise 
time  for  the  blast  wave  and  the  low  frequency  content  found  within  the  acoustic  signals  is  least  attenuated 
over  long  propagation  distances  when  compared  to  some  of  the  predominant  features  initially  identified 
for  discrimination.  From  the  plots  in  Figure  5,  it  can  be  seen  that  at  varying  levels  of  decomposition 
beyond  level  2,  the  energy  distribution  within  the  details  differ  dramatically  prior  to  the  max  peak  pressure 
of  each  blast  with  respect  to  energy  distributions  after  the  main  blast  has  occurred.  Furthermore,  this 
energy  is  not  amplitude  dependant  as  the  baseband  signal  is  filtered  by  the  scaling  function  and  only  the 
high-frequency  noise  components  are  captured  within  the  details  after  decomposition. 

Let  tP  denote  the  sample  time  where  the  maximum  peak  over  pressure  of  a  blast  wave  occurs.  Let 
t0  =  a*fs  designate  an  integration  start  time,  and  tF  =  ft*  fs  denote  integration  stop  time.  We  restrict  the 
values  a,  P  such  that  0  <  a.  p  <  1  for  tP>t0  and  tF  >  t0  such  that  tP ,  t0 ,  tF  e  I .  Define  the  energy 
distribution  within  the  details  just  prior  to  the  max  peak  pressure  as 

1  -If  ' 

°k  =jjYj\D^nl  (3) 

"=‘o 

where  N  =  tP  -  t0  .  Define  the  energy  distribution  at  the  point  immediately  after  the  max  pressure  occurs 
as 

<4> 

m=tP 

where  M  =  tF  -  tP.  Since  Dk ,  Dk  >  0  we  are  able  to  define  analytic  features  xDk  e  SR*  using  the  ratio  of  the 
normalized  energy  distributions  within  the  details. 


x  Dk 


—  log 


10 


Di 

\Dtj 


(5) 
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The  first  3-tuple  of  the  features  space  consists  of  xD5,xD4,xD3  which  we  propose  are  useful  to  characterize 
the  artillery  blasts  of  interest.  A  fourth  feature  is  derived  from  the  wavelet  coefficient  obtained  at  level  5. 
We  integrate  the  magnitude  of  the  area  for  the  coefficients  between  the  start  and  stop  times  defined  for  the 
details 


^SAREA 


=  log 


10 


k=tn 


(6). 


(a)  Oriqnal  HE  signature 


site  5  -  PS  I  Test 


(b)  Approximations  at  level  5  -  DSI  Test 


(cl  Details  at  level  5  -  PS  I  Test 


(d)  Details  at  level  4  -  DSI  Test 


.4Ui.iIu.MjiJ 

(e)  Details  at  level  3  -  DSI  T est ( ; 

Time  (Second) 

Figure  5a  -HE  Signal  and  Wavelet  Decomposition  to  Level  5 


Figure  5b  -Type  A  CB:  Original  Signal  and  Wavelet  Decomposition  to  Level  5 


Time  (Second) 


Figure  5c  -  Type  B  CB:  Original  Signal  and  Wavelet  Decomposition  to  Level  5 
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Figure  6  (a-f)  shows  all  permutations  of  2-D  subspaces  for  the  entire  set  of  data  collected  at  Dugway 
Proving  Ground  that  was  used  to  train  and  initially  benchmark  the  performance  of  a  neural  network 
classifier.  An  ‘x’  is  used  to  indicate  features  derived  from  the  CB  blasts  while  the  corresponding  HE 
feature  points  are  depicted  with  a  Note  that  the  resulting  subspaces  reveal  a  high  degree  of  separability, 
notwithstanding  the  540  -  1660m  variability  in  range  at  which  the  data  was  acquired. 


log  A5  Area  vs  log  D5  Power  Ratio 


log  D5  Power  Ratio  vs  log  D4  Power  Ratio 
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log  D5  Power  Ratio  vs  log  D3  Power  Ratio 
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log  A5  Area  vs  log  D3  Power  Ratio 
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Figure  6  -  2D  Feature  Spaces  derived  from  the  “known”  DPG  data  set  (a)  A5area  vs.  xD3  (b) 
A5area  vs.  xD4  (c)  A5area  vs.  xD4  (d)  xD5  vs.  xD4  (e)  xD5  vs.  xD3  and  (f)  xD4  vs.  xD3 
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3.  NEURAL  NETWORK  CLASSIFIER 

Neural  networks  are  a  powerful  tool  used  to  solve  difficult  classification  (mapping)  problems  with 
complex  decision  regions  that  are  often  required  to  ensure  separability  between  classes  and  a  proven 
ability  to  realize  non-linear  discriminant  functions..  The  use  of  neural  networks  for  classification  is  well 
documented  and  the  requirements  for  training  are  well  known.  The  standard  multilayer  feedforward 
neural  network  shown  in  Figure  7  was  chosen  due  to  its  ability  to  learn  mappings  of  any  complexity, 
provided  that  the  number  of  hidden  layer  neurons  is  sufficient  to  accommodate  the  number  of  separable 
regions  that  are  required  to  solve  a  particular  classification  problem.  In  general,  the  network  contains  Nt 
inputs,  Nh  hidden  layer  neurons  and  n0  output  layer  neurons  with  no  interconnections  within  a  single 
layer.  In  the  three-layer  network  shown,  the  connection  weight  between  the  ith  input  and  jth  hidden  layer 
neuron  will  be  denoted  as  wtj  and  vjk  will  be  used  to  denote  the  connection  weight  between  the  jth 
neuron  in  the  hidden  layer  and  kth  output  layer  neuron.  A  bias  for  the  jth  hidden  layer  neuron  is  denoted 
as  bf  and  the  bias  for  the  k lh  output  neuron  is  bk.  xp  denotes  the  set  of  T ,  A -tuple  feature 
vectors  x„  e '.)?  v  ,  xp  =  [xf  ,x?  J  used  to  discriminate  between  the  different  classes  of  artillery. 

b\  b\ 


Nt  Nh  N0 


Figure  7-  Standard  Multilayer  Feedforward  Neural  Network 

4.1  RESULTS 

Features  were  extracted  from  both  data  sets  using  the  DWT  as  described  in  section  3,  and  the  4-tuple 
feature  vector  xp  =  [x()5,  x,^,  xpm,  a£area\  was  constructed  according  to  equations  (5)  and  (6).  Experiments 
were  then  conducted  to  measure  the  separability  of  the  feature  space  and  benchmark  the  performance  of 
the  proposed  neural  network  classifier.  The  architecture  shown  in  Figure  7,  consisting  of  a  single  hidden 
layer  neuron  was  trained  using  22  randomly  selected  vectors  from  a  total  of  46 1  signatures  made  available 
from  the  data  collection  exercise  at  DPG.  The  training  set  comprised  of  1 1  samples  of  HE  data  and  1 1 
samples  of  Type  A  CB  rounds  collected  at  each  of  the  3  sensor  sites.  The  neural  network  was  trained  with 
sample  features  representing  a  HE  round  would  result  as  zero,  while  when  features  from  a  sample  CB 
round  applied  would  output  a  one  based  on  the  properties  of  our  neural  network.  The  network  was  trained 
with  the  total  error  of  the  system  being  less  than  5e-3  and  a  learning  rate  of  0.1  was  utilized.  When  tested 
against  the  remaining  439  signatures,  the  networked  correctly  classified  100%  of  the  remaining  225  CB 
blasts  and  misclassified  a  total  of  4  out  214  (or  98.1%)  of  the  remaining  HE  rounds,  for  an  overall 
performance  of  99.1%.  The  rounds  were  classified  based  on  the  neural  network  generated  weights  as  CB 
if  the  output  value  was  greater  than  lA  and  HE  otherwise. 
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In  a  subsequent  experiment,  a  neural  network  consisting  of  4  hidden  layer  neurons  was  then  trained  using 
the  entire  dataset  from  the  DPG  test;  a  total  of  236  Type  A  CB  rounds  and  225  HE  signatures.  The 
network  was  then  tested  against  the  data  collected  by  NCPA  at  YPG.  The  network  correctly  classified  165 
of  the  total  of  166  CB  rounds  and  misclassified  a  total  of  6  out  of  57  HE  rounds,  for  a  rate  of  99.4%  and 
89.5%  respectively.  Neural  networks  consisting  of  up  to  3  hidden  layer  neurons  could  not  match  or 
surpass  the  classification  performance  provided  with  the  architecture  comprised  of  4  hidden  layer  neurons. 
A  third  experiment  was  undertaking  utilizing  the  “so-called”  blind  data  from  DPG  for  which  truth  data 
was  provided  by  the  test  facility.  The  data  was  tested  using  a  neural  network  trained  with  4  hidden  layer 
neurons  and  against  all  the  initial  DPG  data  that  was  considered  to  be  “not  blind”.  The  results  obtained  for 
these  blasts  provided  correct  classification  results  of  98.3%  of  all  the  CB  and  95.7%  of  all  the  HE, 
providing  further  confidence  in  the  neural  network  and  the  results  of  the  classification.  However,  the 
same  results  were  obtained  when  training  a  neural  network  comprised  of  5  or  more  hidden  layer  neurons. 
Once  again,  the  rounds  were  classified,  as  CB  when  the  value  produced  was  greater  than  V2  and  HE 
otherwise.  A  summary  of  the  neural  network  performance  is  shown  in  Table  1. 


Table  1  -  Summary  of  Classification  performance  results 


Experiment 

# 

Training  Data 

Test  Data 

Classification 

Percentage 

1 

1 1  CB  (DPG) 

225  CB  (DPG) 

225  CB  /  0  HE 

100% 

1 1  HE  (DPG) 

214  HE  (DPG) 

210  HE/ 4  CB 

98.1  % 

2 

236  CB  (DPG) 

166  CB  (YPG) 

165  CB  /  1  HE 

99.4  % 

225  HE  (DPG) 

57  HE  (YPG) 

51  HE/6CB 

89.5  % 

3 

236  CB  (Blind) 

230  CB  (Blind) 

226  CB  /  4  HE 

98.3  % 

225  HE  (Blind) 

184  HE  (Blind) 

174  HE/8  CB 

95.7  % 

4.2  REAL  TIME  PERFORMANCE 

The  algorithm  described  in  this  paper  was  coded  into  a  real-time  environment  using  C++  and  incoiporated 
into  an  application  for  an  acoustic  based  CBRN  discrimination  system  into  an  unattended  ground  sensor 
(UGS)  upon  the  PAWSS  platform.  The  algorithm  was  embedded  into  PAWSS  as  part  of  a  LOE  to 
provide  processed  acoustic  event  data  as  added  situational  awareness  to  the  user  and  the  decision  aide  tool. 
A  result  is  provided  to  the  user  via  a  graphical  user  interface  (GUI)  that  indicates  whether  a  burst  event 
was  a  HE  or  potentially  CB  in  nature.  The  sensor  layout  is  shown  in  Figure  8.  Microphones  were 
emplaced  at  positions  S1-S4  at  sampling  rates  of  32KHz,  32kz,  16kHz  and  6.25  kHz  respectively.  S3  was 
emplaced  for  the  purpose  of  recording  the  muzzle  blast,  S2  was  used  to  record  impact  events  and  S4  was 
used  to  test  the  real  time  implementation  of  the  algorithm. 
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The  results  obtained  for  38  airburst  events  resulted  in  100%  correct  classification  of  the  CB  and  HE  single 
volley  events  at  the  ranges  shown  in  Figure  8.  The  detonation  ranges,  airburst  heights,  ever-changing 
meteorological  conditions,  and  artillery  stock  varying  greatly  during  the  3-day  test  period  provided  a  host 
of  challenges  for  the  algorithm  classifying  single  volleys.  The  PAWSS  LOE  provided  an  opportunity  to 
benchmark  the  algorithm  against  dual  shot  volleys  for  the  first  time.  A  total  of  34  total  airburst 
detonations  were  used  to  test  against  the  dual  shot  volleys.  The  dual  shot  volleys  were  designed  to 
simulate  a  large  scale  attack  by  simultaneously  firing  both  howitzers  followed  immediately  by  a  third 
singular  round  in  a  series  of  combinations.  All  possible  arrangements  for  the  3  rounds  of  the  two  artillery 
types  were  employed  (i.e.  HE,  HE,  CB  or  HE,  CB,  HE  or  HE,  CB,  CB,  or  CB,  CB,  HE)  to  both  mimic 
field  conditions  and  mask  airburst  events.  The  algorithm  classified  80%  of  all  the  HE  rounds  and  83.3% 
of  all  the  CB  events  (cf.  Table  2). 


Table  2  -  Summary  of  Classification  performance  results 


Experiment  #  4 

Training  Data 

Test  Data 

Classification 

Percentage 

PAWSS  LOE  Single  Shot 

1 1  CB  (DPG) 

24  CB  (DPG) 

24  CB  /  0  HE 

100% 

1 1  HE  (DPG) 

14  HE  (DPG) 

14HE/0  CB 

100% 

PAWSS  LOE  Dual  Shot 

1 1  CB  (DPG) 

24  CB  (DPG) 

24  CB  /  4  HE 

83.3% 

1 1  HE  (DPG) 

10  HE  (DPG) 

8  HE / 2  CB 

80% 

5.  CONCLUSION 

In  this  paper,  feature  extraction  methods  based  on  the  discrete  wavelet  transform  and  multiresolution 
analysis  facilitate  the  development  of  a  robust  classification  algorithm  that  affords  reliable  discrimination 
between  conventional  and  simulated  chemical/biological  artillery  rounds  via  acoustic  signals  produced 
during  detonation.  Amplitude  dependant  features  provide  the  basis  for  previously  developed 
discrimination  algorithms,  which  vary  dramatically  due  to  signal  attenuation  and  distortion  resulting  from 
various  propagation  effects.  Measurements  from  multiple  sensors  can  be  used  to  obtain  range  information 
to  a  blast  under  these  circumstances,  which  in  turn  maintains  the  potential  for  a  high  level  of  correct 
classification.  In  general,  it  was  noted  that  when  the  distances  from  the  detonation  point  to  disparate 
sensor  exceeds  1km;  the  probability  of  correct  classification  degrades  dramatically,  even  though  visual 
cues  that  can  be  used  to  exploit  class  specific  attributes  are  still  noticeable  within  the  signature  data. 

The  non-stationary,  transient  and  often  oscillatory  nature  of  artillery  blast  signals  may  be  characterized 
efficiently  with  wavelet  bases  that  successfully  capture  the  time-frequency  distribution  of  such  signal 
components.  Wavelets  are  better  suited  for  analyzing  transient  signals,  as  they  are  well  localized  in  time 
and  in  frequency,  and  are  able  to  take  into  account  the  scale  of  discrete  signal  components.  To  this  end, 
the  wavelet  transform  has  been  shown  to  provide  a  scalable  time-frequency  depiction  of  blast  signatures 
and  has  revealed  details  that  are  not  readily  detectable  using  conventional  signal  processing  techniques. 
Utilizing  the  idea  of  multiresolution  analysis  provided  by  the  discrete  wavelet  transform,  we  were  able  to 
break  the  original  signature  into  subband  components  thereby  removing  the  higher  frequency  noise 
features  and  creating  two  sets  of  coefficients  at  varying  levels  of  decomposition.  These  coefficients  are 
obtained  each  time  the  signal  is  passed  through  a  lowpass  and  highpass  filter  bank  whose  impulse 
response  is  derived  from  Daubechies  db5  wavelet.  Distinct  features  are  obtained  through  the  process  of 
isolating  the  details  of  the  higher  oscillatory  components  of  the  signature  (i.e.  characteristic  whining  of  the 
shrapnel  for  HE  rounds  and  weak  under  pressure  for  CB  rounds).  The  ratio  of  the  energy  contained  within 
the  details  at  varying  levels  of  decomposition  is  sufficient  in  discriminating  between  the  vast  majorities  of 
artillery  blasts. 

The  discrete  wavelet  transform  has  been  successfully  employed  to  extract  unique,  disjoint  feature  sets  that 
remain  consistent  for  a  given  class  of  round  and  do  not  degrade  dramatically  with  long-range  propagation. 
Classification  results  using  a  neural  network  trained  on  such  attributes  show  consistently  reliable 


RTO-MP-SET-107 


8  - 11 


UNCLASSIFIED/UNLIMITED 


UNCLASSIFIED/UNLIMITED 


Reliable  Classification  of  High  Explosive  and 
Chemical/Biological  Artillery  Using  Acoustic  Sensors 


ORGANIZATION 


discrimination  of  CB  artillery  rounds  that  surpasses  98%,  even  when  the  network  is  trained  on  data 
acquired  during  different  field  test  experiments.  In  addition,  the  algorithms  developed  herein  allow  for  the 
discrimination  of  base  ejection  type  rounds,  which  were  not  part  of  the  original  training  set  used  to  train 
the  neural  network. 

The  algorithm  developed  herein  were  implemented  at  the  PAWSS  LOE  there -by  validating  the  distinct 
feature  set  extracted  from  the  airburst  events  discriminating  CB  and  HE  detonated  events.  The  feature  set 
remains  consistent  for  a  given  class  of  rounds  (independent  of  payloads)  and  different  field  test 
experiments.  Test  results  show  that  these  features  do  not  degrade  dramatically  as  long-range  propagation 
is  increased.  The  real  time  implementation  was  tested  at  increased  distances  to  the  impact  zone;  upwards 
of  3000m  with  burst  heights  of  upwards  of  1000m.  The  system  was  put  through  several  extreme 
conditions  common  to  hostile  environments  and  performed  commendably  against  high  wind  gusts  in 
excess  of  30m.p.h.,  rain,  rugged  mountainous  terrain  and  multiple  shot  volleys.  The  results  of  the  PAWSS 
LOE  test  were  100%  correct  classification  for  single  volley  rounds,  which  the  system  was  trained  to  solve. 
The  algorithm  trained  only  against  single  volleys  for  the  first  time  was  exposed  and  tested  against  multiple 
volley  rounds.  The  system  correctly  identified  83%  of  the  events  providing  impressive  real-time  results 
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Purpose  /  Motivation 
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Chemical  and  Biological  Weapon 

Threats  and  Needs 


•  Determining  if  explosive  event  contains  High 
Explosive  material  or  Chemical/Biological  agent 
on  the  battlefield. 

•  Providing  field  commanders  with  greater  response 
time  using  a  stand  alone  acoustic  sensor. 

•  Giving  greater  situational  awareness  to  threatened 
soldiers. 
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Data  Collection 
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Acoustic  Signature  Data 
Collection  of  Blast  Events 


•  Data  Collection  1 . 

Conducted  by  National  Center  of  Physical  Acoustics  (NCPA)  in 
cooperation  with  ARDEC. 

-  39,  rounds  fired. 

3  categories  of  rounds  were  used,  HE,  Type  A  CB,  and  Type  B. 

•  Data  Collection  2. 

Conducted  by  DPG  Team  and  U.S.  Army  Edgewood  Chemical 
Biological  Center  (ECBC) . 

-  265,  rounds  fired. 

-  2  categories  of  rounds  were  used,  HE  and  Type  A  CB. 
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Data  Collection  #1 


A  C  0  E  Data  C  ollection 


Detonation  Impacts  ite 
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Typical  Blast  of  HE  Round 


•  High  frequency  precursors 
to  the  main  blast. 

Generated  by  Supersonic 
Shrapnel  Elements. 

•  Large  Amplitude  of  Main 
Blast. 

•  Large  under  pressure 
element . 

-  Generated  by  large 
comparable  weight  of 
explosives  rapidly  burning. 

•  Short  Duration  Signatures. 


HE  Round 
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Typical  Blast  of  Type  A  CB  Round 

•  Small  amplitude 
associated  with  main  blast. 

-  The  explosive  material  is 
minimal  compared  to  the 
comparable  HE  round  type. 

•  Elongated  bum  time 
following  main  blast. 

-  The  deliberately  slow  to 
properly  release  the 
compounds. 

•  Weak  under  pressure. 
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Typical  Blast  c 


•  Short  Duration  Pulse. 

-  Resulting  from  base 
ejection  rounds. 

•  Weak  Under  Pressure. 

-  Small  amount  of 
Explosives. 

•  Slow  Bum  Time. 

-  Elongated  to  properly 
discharge  contents  of 
the  round. 
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Signal  Processing 
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Discrete  Wavelet  Transform 

(DWT) 

•  Derived  from  subband 
filters  and  multiresolution 
decomposition. 

-  Coarser  Approximation. 

-  Removing  high  frequency 
detail  at  each  level  of 
decomposition. 

•  Acts  like  a  multiresolution 
transform. 

-  Maps  low  frequency 
approximation  in  coarse 
subspace  high  frequency 
elements  in  a  separate 
subspace. 
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Daubechies  Wavelet,  n  =  5 


•  Representation  of  the  scaling  and  translation  function  of 
db5. 

Scaling  function  resembles  blast  signature  of  the  HE  and  CB 
rounds. 

Provides  the  ability  to  approximate  signal  with  the  characteristic 
wavelet. 
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Multiresolutional  Analysis 


•  Using  a  series  of  successive  high  pass  and  low  pass  filters  to 
create  a  set  of  subspaces. 

-  High  pass  filter  obtains  the  details  of  the  signatures  while  the  low 
pass  filter  obtains  a  coarse  approximation  of  the  signal. 

•  The  resulting  banks  of  dyadic  multirate  filters  separate  the 
frequency  components  into  different  subbands. 

-  Each  pass  through  gives  you  resolution  of  factor  2. 
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Effects 


Wavelet  decomposition  to 
level  5  of  three  varying 
blast  types  from  varying 
ranges. 
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Wavelet  Extracted  Features 


Comprised  of  primitives 
derived  from  the  normalized 
energy  distributions  within  the 
details  at  level  5,  4,  and  3  of  the 
wavelet  decomposition. 

Distribution  of  blast  type  differ 
greatly  when  takeiymoi^c^he 
max  pressure, 
with  respect  to  distribution  after 
the  max  blast, 


Resulting  Ratio 


( 


XDk  =1°gl 


A5  area  is  a  feature  derived 
from  wavelet  coefficients  at 
level  5. 

Integrating  the  magnitude  of  the 
area  for  the  coefficients 
between  the  start  and  stop 
times. 
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Extracted  Features  Using  DWT 
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4- 


This  energy  ratio  leads  to 
the  discover  of  4  features 
with  A5  area  that  are  not 
amplitude  dependent. 


Our  n-tuple  feature  space 
thus  becomes  a  4-tuple 
space,  ,  ti 


be  applied  for 
classification. 
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2-D  Feature 


log  A5  Area  vs  log  D5  Power  Ratio 
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Space  Realization 


log  D5  Power  Ratio  vs  log  D4  Power  Ratio 
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Airburst  Discrimination 
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Neural  Network 


•  Realize  non-linear 
discriminant  functions  and 
complex  decision  regions 
to  ensure  separability 
between  classes. 

•  Standard  Multilayer 
Feedforward  Neural 
Network. 

•  Number  of  hidden  layer 
neurons  depend  on 
complexity  of  required 
mapping. 
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Results  of  Training  Neural 


Network  to  DSI  Data 


Feature  Space  created 
using  DWT. 

4-tuple  feature  vector. 


p  _  P  P  P  AP 
x  ~  XD5’XD4’XD3’A5AREA 


22  randomly  selected 
vectors  from  46 1 
signatures. 

Trained  Neural  Network 
to  trained  output  data  of  0. 

Single  hidden  layer  neuron. 

Total  error  in  equation  after 
training  is  less  then  5e-3. 

-  Learning  rate  of  0. 1 . 
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Results  of  HE/CB  Discrimination 


•  Experiment  1. 

Applying  a  neural  network 
with  the  weights  in  the  table 
1  to  DPG  data,  99.1% 
Correct  Classification. 

•  Experiment  2. 

-  A  neural  network 

containing  4  hidden  layer 
neurons  trained  using  entire 
DPG  dataset  tested  against 
NCPA  dataset,  96.9% 
Correct  Classification. 


Experiment  # 

Training 

Data 

Test  Data 

Classification 

Percentage 

1 

11  CB 

(DSI) 

225  CB 

(DSI) 

225  CB  /0  HE 

100% 

11  HE 

(DSI) 

214  HE 

(DSI) 

210  HE  /  4  CB 

98.10% 

2 

236  CB 

(DSI) 

166  CB 

(YPG) 

165  CB /I  HE 

99.40% 

225  HE 

(DSI) 

57  HE 

(YPG) 

51  HE /6  CB 

89.50% 

% 

Wi? 

W/4 

vn 

11.6967 

0.5343 

-0.4958 

-2.4991 

-13.4966 

4.6377 

1.2455 

3.5569 

5.3068 

13.3761 

4.7023 

0.9875 

7.3951 

8.902 

-15.3761 

-5.2246 

1.481 

2.6982 

4.1203 

-19.6513 

-2.8169 

1.4847 

-18.9732 

-23.6088 

-14.286 
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Blind  Results  of  HE/CB  discrimination 


•  Experiments. 

-  Utilizing  the  neural  network 
containing  4  hidden  layers 
neurons  trained  against  the 
entire  “known”  DPG  data 
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wn 

wj4 
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7.3951 
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-15.3761 

-5.2246 

1.481 

2.6982 

4.1203 

-19.6513 

-2.8169 

1.4847 

-18.9732 

-23.6088 

-14.286 

set  was  then  tested  against 
the  “blind  data”  the  results 
once  compared  with  the 
truth  resulted  in  98.3%  and 
95.7%  reliable 
classification. 
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Real-Time  Implementation 
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Experiment  4  Real  Time  Implementation 


Portable  Area  Warning  Surveillance 
System  (PAWSS). 

lyr  Limited  Objective  Experiment  (LOE). 

-  Focused  on  the  utility  of  cascading 
detection  methodologies. 

Combines  Stand-off  CBRN  systems  to 
address  both  force/installation  protection. 

LOE  Outcomes. 

-  Operable  Products  leading  to  fully 
designed  products  that  are  sustainable. 

Demonstration  of  capabilities  within 
simulated  battlefield  environments  of 
layered  wide  area  cascading  detection. 
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PA  WSS  LOE  Test  Layout 
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PA  WSS  LOE  Results 

•  June  1 9th-28th  Portable  Area  Warning  Surveillance  System 
(PAWSS)  Limited  Objective  Experiment  (LOE). 

•  Implemented  real  time  version  of  CBRN  Discrimination  at 
PAWSS  LOE  conducted  by  ECBC. 

•  1 00%  single  volley  discrimination,  never  tested  against 
dual  volley,  still  83%,  also  all  event  starts  were  detected 
for  100%. 

•  Assist  in  transition  and  support  of  acoustic  element 
CBRNEWS  ATD  extending  LOE  efforts. 


Event  Type 

#  of  Events 

Discriminated  Correctly 

Single  Round 

38 

38/38;  100% 

Dual  Round 

34 

28/34;  83% 
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Real  Time  Performance 


•  During  June  21st  and  June  22nd,  2005  a  proof  of  concept  test  was  conducted  for  the 
acoustic  CBRN  discrimination  algorithm. 

-  PAWSS  Test  Site,  DPG. 

•  Acoustic  System  from  Impact  Zone. 

-  A  C++,  real  time  algorithm  was  tested  at  DPG  as  part  of  the  acoustic  portion  of  PAWSS  LOE 
conducted  by  JPM  for  NBC  Contamination  Avoidance  at  ECBC. 

-  A  total  of  72  HE/CB  rounds  were  detonated. 

•  A  howitzer  fired  24  HE,  and  48  CB  rounds. 

•  Single  Round  Volley  Results. 

38  Airburst  Detonation  (14  HE,  24  CB),  Correct  Classification. 

•  Multiple  Round  Volley. 

-  CBRN  Algorithm  Never  Benchmarked  in  Lab  vs.  Multiple  Rounds. 

•  2  Rounds  simultaneously  fired  followed  by  a  3rd  round  fired  soon  as  possible. 

-  34  Airburst  Detonation  (10  HE,  24  CB). 

•  17  events,  each  event  consisted  of  2  detonations. 

-  83%  Overall  Correct  Discrimination  of  HE/CB. 

•  discrimination  on  all  HE  rounds. 

•  acoustic  detection  of  all  events. 

•  28  correctly  discriminated  from  34  detonations. 

•  Shortcomings  occur  within  the  data  acquisition  process,  limited  by  processing  window  size. 
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Conclusion 

•  Features  extracted  facilitate  robust  classification. 

-  Reliable  discrimination  of  CB  rounds,  or  greater  of  single  volley  events. 

•  The  features  this  algorithm  is  based  on  go  beyond  previous  amplitude  dependent 
features. 

-  Degradation  due  to  signal  attenuation  and  distortion  is  nullified  and  exceeds  3km  in  range  propagation. 

•  Scalable  time  frequency  representation  uncovered  non-readily  detectable  features. 

-  Subband  components  remove  higher  frequency  noise  features. 

-  Isolating  the  details  of  higher  oscillatory  components. 

•  Real  time  verification  at  PAWSS  LOE  of  CBRN  Discrimination  Program 
Implemented  in  C++. 

-  Single  volley  round  discrimination  in  real  time  for  all  variants  was 

-  Dual  volley  round  discrimination  in  real  time  for  all  variants  was  ,  and  detected  an  event  of 
the  time. 

•  Wavelets  can  be  possibly  used  to  discriminate  varying  types  of  artillery  projectile 
launches  from  impacts  independent  of  range. 

-  Utilizing  wavelets  and  other  signal  processing  techniques  to  perform  a  similar  task  as  described  within 
with  refinement  for  the  problem. 

•  Future  Considerations. 

-  Networking  of  sensors  can  provide  TDOA  abilities  to  further  localize  a  threat. 

-  Transitioning  algorithm  to  CBNEWS  ATD. 


3/27/2008 


SET-1 07/RYS-020  Session  2,  #8 


30 


QUESTIONS...? 
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BACKUP 
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•  Efficiently  represent  non¬ 
stationary,  transient,  and 
oscillatory  signals. 

•  Desirable  localization 
properties  in  both  time 
and  frequency  that  has 
appropriate  decay  in  both 
properties. 

•  Provide  a  scalable  time- 
frequency  representation 
of  artillery  blast  signature 


3/27/2008 


SET- 


Scale 


Wavelets 


Time 


-020  Session  2,  #8 


34 


